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ABSTRACT 



We examine the velocity structure in the gas associated with HI in the 
damped Lya absorption system at redshift z = 1.7764 towards the QSO 
1331 + 170 using Arecibo HI 21cm data, optical spectra from Keck HIRES and 
ESO VLT-UVES, and a previously published HST STIS ultraviolet spectrum. 
From the optical data we find at least two, and possibly three, components 
showing CI lines. One of these has very narrow lines with Doppler parameter 
b = 0.55km s" 1 , corresponding to a kinetic temperature of 220K if the broadening 
is thermal and with a 2-a upp er limit of 480K. We re-examine the H 2 analysis 



undertaken by 



Cui et al. 



( 120051 ) using the neutral carbon velocity structure, and 
find a model which is, unlike theirs, consistent with a mixture of collisional and 
background radiation excitation of the observed H 2 rotational levels. 

Using Voigt profile fits to absorption lines from a range of singly ionized 
heavy elements we find eight components covering a velocity range of ~ 110 
km s -1 , with a further outlier over 120 km s -1 away from the nearest in the main 
group. The HI structure is expected to follow some combination of the singly 
ionized and neutral gas, but the 21cm absorption profile is considerably different. 
We suggest, as have others, that this may be because the different extent and 
brightness distributions of the radio and optical background sources mean that 
the sightlines are not the same, and so the spin temperature derived by comparing 
the Lya and 21cm line strengths has little physical meaning. The neutral and 
singly ionized heavy element line profiles also show significant differences, and so 
the dominant components in each appear to be physically distinct. Attempts to 
use the range of atomic masses to separate thermal and turbulent components of 
their Doppler widths were not generally successful, since there are several blended 
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components and the useful mass range (about a factor of two) is not very large. 

The velocity structure in all ionization stages up to +3, apart from the neutral 
heavy elements, is sufficiently complex that it is difficult to separate out the 
corresponding velocity components for different ionization levels and determine 
their column densities. 



Subject headings: Galaxies: ISM — quasars: absorption lines — quasars: individual 
(1331+170) 
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Introduction 



It has long been known that there is significant temperature structure in the interste llar 
mediu m in the Galaxy. Theoretical models, as originally proposed by e.g. 



Field et al 



( 119691 ) and 



McKee &: Ostrikerl ( 119771 ) and refined by several authors subsequently, predict 
a number of phases in pressure equilibrium. These are a cold neutral medium (CNM), 
characterized by temperatures T of up to a few 100 K; warm neutral and warm ionized 
media (WNM and WIM) with 1000 < T < 10000 K, and a hot ionized medium (HIM) with 
T > 10 6 K. 

Within the Galaxy, observations of the 21 cm line have led to the detection of hyperfine 
spin temperatures ~ 3000 - 8000 K, characteristic of the warm neutral medium (Carilli, 



Dwarakanath & Goss, 1998, Kanekar et a. 



2003). Observations of both the CNM and 



WNM are described by 
the cold medium up to ~ 200 K, and upper 
~ 20 000 K. Contrary to the predictions of 



Heiles fc Trolandl (120051 ). who give a range of spin temperatures for 



imit kinetic temperatures for the WNM of up to 



McKee fc Ostrikerl (119771 ) a significant fraction 



of the WNM is at temperatures where the gas is thermally unstable; i.e. T w 500 to 1500 
K . On the other hand, direct observational confir mation of the temperature s predicted 



Redfield k Linskvi (120041 ) . They used 



for a thermally stable WNM have been provided by 
high resolution spectra from the Hubble Space Telescope and compared line profiles from 
DI with those from heavy elements including Mgll and Fell to obtain local interstellar 
medium temperatures ~ 6500 ± 1500 K. They also noted that the line broadening is not 
always purely thermal, and that there is a mean bulk flow component within regions of a 
little over ~ 2 km s . 

For most heavy element quasar absorption systems the gas is ionized and the kinetic 
temperature estimates range from ~ 10 4 K and up to ~ 10 6 K in some cases where O VI 
is found. For those with high HI column densities, the damped Lya systems (DLA), the 
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gas is clearly neutral, with the heavy elements neutral or singly ionized. However, the 



ubiquitous presen ce of CIV and Si IV ( I Wolfe &: Prochaska 



of O VI and N V flFox et al. 



2000l ) and the recent detection 



20071 ) absorption indicates that these high redshift gas layers 
are also multi-phase media. The HI 21 cm spin temperatures provide a means of estimating 
the thermal conditions in the gas, and generally spin temperatures are measured to be 
T s > 10 3 K (see e.g. the compilation by Curran et al., 2005). On the other hand these 
estimates assume that the HI column density, iV(Hl), inferred from the Lyct absorption 
can be used to compute T s from the 21 cm optical depth. While this is a straightforward 
procedure which has been used for the past 30 years (see Wolfe & Davis, 1979), it is 
complicated by the disparity between the physical sizes of the optical and radio beams 
subtended by the background quasar at the absorption redshift. Whereas the optical beams 
are typically less than one light year across, VLBI experiments reveal that the radio beams 
are typically a few hundred pc at the low frequencies of the radio absorption lines (e.g. 
Briggs et al., 1989; Polatidis et al.,1995; Thakkar et al., 1995). Moreover, the appearance of 
shifts in interferometer phase between the 21cm line and continuum in VLBI experiments 
demonstrates a shift in the brightness centroid of the background radio source, which can 
only be reproduced by a non-uniform distribution of 21 cm opacity toward radio sources 
with scale-lengths of over 200 pc (e.g. Wolfe et al., 1976). For these reasons, the values of 
T s inferred by these techniques should be regarded as conservative upper limits. 

Temperatures may also be inferred from excitation conditions in other ions and from 
molecula r hydrogen. For example, u sing m easurements of the fine structure levels of Sill 



and C II iHowk. Wolfe fc Prochaska! (120051 ) found that the excitation conditions constrain 
the temperature to T < 954 K in one high redshift DLA. In those DLAs with H 2 absorption, 
analyses of the relative populations of the lower rotational states give kinetic temperatures, 
which are generally less than ~ 200 K (e.g Cui et al., 2005; Ledoux, Petitjean & Srianand, 
2006). 
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The relationship between excitation temperatures T ex and kinetic temperatures, T, 
depends on the balance between collisional and radiative processes. In DLAs exhibiting 21 
cm absorption the hyperfine excitation temperature, i.e. the spin temperature T s , likely 
equals the gas kinetic T because the hyperfine levels are populated mainly by collisions 
owing to the moderately high densities for the gas, and the long radiative lifetime of the 
upper (.F=l) hyperfine state. Furthermore, even in regions with low density, ambient Lya 
radiation can couple T s to T via the Field- Wouthuysen effect (e.g. Field, 1959). Within 
the Galaxy, there is good agreement between T and T ex determined from the relative 
populations of the J=0 and J—l rotational states of H 2 , an agreement that holds for 
column densities logiV(H 2 ) <; 16 (cm -2 ). However, excitation temperatures deduced from 
the relative populations of levels with J > 2 depart from T since spontaneous photon 
emission rates exceed collisional de-excitation rates due to the short radiative lifetimes. 
Rather the values of T ex inferred for these states is determined by the ambient radiation 
fields and other processes responsible such as H 2 formation on grain surfaces. 

In this paper we obtain upper limits to the gas kinetic temperature in neutral regions 
using the CI absorption lines arising in a DLA at redshift z = 1.77642 to wards the quasa r 



Q1331 + 170. Because this DLA has been detected in 21cm absorption (j Wolfe fc Davis 



19791 ). we can compare this temperature with the spin temperatur e of T s > 1000 K. Also 



20051 ) has been obtained for 



the H 2 excitation temperature T ex = T H2 = 150 K ( ICui et al. 
this DLA. We shall argue that the true kinetic temperature of some of the gas is at most a 
few 100 K, and that the value of the spin temperature is likely to be misleading as a result 
of non-uniform coverage of the background radio source. We also re-examine the physical 
processes responsible for thermal equilibrium in this DLA and remark on a dilemma 
not previously d iscussed. Namely, for a plausible range of volume densities, background 



radiation alone ( IHaardt fc Madau 



200JJ) will drive T ex » T. For this reason, the reported 



H 2 level populations in this DLA are difficult to understand. However, adopting the 
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velocity structure inferred from the C I and applying it to the H2 does allow us to produce 
a physically consistent picture. 



We also examine the singly ionized heavy elements, and two more highly ionized 
species, to look for correspon ding structure. This has be en d one for several objec t s as p art 



of large DLA surveys e.g. by 



Prochaska fe Wolfe! (|1999h and 



Wolfe & Prochaskal (|2000|), 



with the general conclusion that the singly ionized species structure is similar to that of 
A1III, while the component structure of e.g. CIV and singly ionized species bear little 
relation to each other. We find generally similar results, and show that, for DLA1331 + 170, 
the fitting of complex velocity structures and associating components across ionization levels 
is not always straightforward. Consequently there wi ll be uncertainties in the ionization 



equilibrium models of DLAs (jHowk fc Sembach 



19991 ). Also, the Mil/ Mill indicator has 



been used for the analysis of sub-DLA systems (e.g. Dessauges-Zavadsky et al., 2003). 



Observations 



2.1. Radio 



The hitherto unpublished radio spectrum was acquired with the Arecibo 300 m 
telescope of the National Astronomy and Ionosphere Center near Arecibo, Puerto Rico 
by one of us (AMW) during April, 1991. Because of improvements in receiver technology, 
the signal-to- noise rat ios (S/N) of th e data were significantly improved over the Arecibo 



discovery spectrum of 



Wolfe fc Davisl ( 119791 ). The S/N for the data varies between 200 and 



400 per 0.7 kms 1 pixel, and the resolution after Hanning smoothing is about 1.4 kms . 
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2.2. Optical 



The Keck HIRES f|Vogt et alJll994h spectrum of Q1331 + 170 is the one described by 



Prochaska Sz Wolf el ( 1997 ). It covers the spectral range from 4220 - 6640 A at a resolution 



(full width half maximum) of 6.25 km s 1 , with a number of gaps at wavelengths > 5515 A. 
The S/N~ 100 per 2 km s -1 pixel at 5400 A, decreasing to ~ 60 in the 4300 - 4600 A region. 



The VLT UVES flDekker et al. 



20001 ) spectrum is based on data available in the ESO 
archive in July, 2006. This was obtained with two of the standard echelle settings: on 
2002-10-04, for ESO program 67.A-0022(A), 3x4500s exposures with blue arm central 
wavelength 346nm and red arm central wavelength 860nm, and on 2003-03-12, for ESO 
program 68.A-0170(A), 2 x 3600s eposures with the blue 390/red 564 setting. 2a;2 on-chip 
binning was used in both cases. See the VLT UVES handbook]] for more details. The data 
were extracted and combined using the UVES_popler package^. The combined spectrum 
covers the range from 3050 A to 10080 A, with gaps at 4517 - 4623, 5597 - 5679, and some 
more at wavelengths > 6650 A. The resolution is 7 km s _1 , and the S/N~ 40 per 2.5 km s^ 1 
pixel at 5400 A. At 4300 A S/N~ 30, but it is better at shorter wavelengths, ~ 35 at 3600 A. 



Voigt profi 



es convolved with the instrument profile were fitted to the data using the 



VPFIT packagd_|, version 9.5. In versions 9 and later the model profiles are generated 
on a finer grid than the data, and in the application here at least 9 sample points were 
required across the intrinsic FWHM of the narrowest line so that the intrinsic line profile 
was adequately modelled. The convolution with the instrumental profile is done before 
re-sampling back to the original pixel scale for direct comparison with the original data, 



1 http:/ /www. eso.org/sci/facilities/paranal/instruments/uves/doc/index.html 
2 http:/ /astronomy.swin. edu.au/ ~mmurphy/UVES_popler.html 
3 http: / /www. ast.cam.ac.uk/ ~rfc/ vpfit.html 
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so there should be no numerical problems associated with possible undersampling of the 
absorption profile which might have arisen with previous versions of the program. 



3. HI spin temperature 



While there appears to be some structure in the 21 cm line profile (see Figs [T] and 
E]), it is at a marginal level. The overall profile is well fitted by a single Gaussian with 
z = 1.7764610 ± 0.0000039, b = 15.7 ± 0.7km s^ 1 . Integrating over the optical depth of the 



entire 21 cm feature gives f r u dV = 0.93 ± 0.03 km 



The HI column density inferred from a Voigt profile fit to the damped Lya line 
for the whole complex fitted as a single component (with an additional component in 
the long wavelength base to allow for the kno wn sub-DLA system at z = 1.78636) is 



logiV(Hl) = 21.17 (cm 2 ), in agreement with iProchaska fc Wolfd (119991 ) who report 



logiV(Hl) = 21.176 ± 0.041. The Lya and 21cm profiles are shown on the same velocity 
scale in Fig[TJ The formal error in our coloumn density value is < 0.01, but systematic 
errors are likely to predominate. Estimates using different continua and blending from other 
systems gave a range 21.10 < logiV(HT)< 21.24, so we adopt 0.07 as an error estimate. 
The redshift is z = 1.77674 ± 0.00009, which is consistent with the 21cm redshift. If we 
constrain the Lya redshift and Doppler parameter to be the same as those for the 21 cm 
line, then logiV(Hl)= 21.17. 

With this value for the HI column density, the spin temperature T s = 870 ± 160 K, 



assuming com 



Wolfe k Davis 



ete c overage of the absorber. This is consistent with the results of 



19791 ) . who give T s = 980 K, with a lower limit of 770 K. If we take the 



covering factor estimate / = 0.72 provided by lKanekar et al. 
temperature is T s = 1360 K. 



)9) 



(2009), then the corresponding 
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4. Neutral Heavy Elements 

We first investigate transitions arising from the neutral state of various elements, since 
these are expected to arise in the CNM phase of the gas. The main lines available in the 
spectral range covered are from CI and Mgl. 



CI in this system has been studied by 



Songaila et al. 



jmj) , using Keck HIRES 



spectra of the ground state and fine structure lines in the multiplets at 1656 and 1560 A to 
determine excitation temperatures of T ex = 10.4 ± 0.5 for the component at z = 1.77638, 
and T ex = 7.4 ± 0.8 at z = 1.77654. The UVES spectra extend the coverage so that lines at 
shorter wavelengths, specifically 1328, 1280 and 1277 A, may also be used to constrain the 
component parameters. Some of these are blended with lines of other elements at different 
redshifts, and all lie within the Lya forest absorption region. We have carefully fit the 
CI lines and included possible blends in the fitting procedure, taking advantage of other 
transitions from these blended lines to help constrain their parameters. Details for each 
wavelength region used in the fitting procedure are as follows: 



• CI 1656 falls in a gap in the VLT UVES coverage, and so the fit relies on Keck 
HIRES data alone. There is blending with weak CIV 1550 at z — 1.96616 in the blue 
wing. C IV 1548 at this redshift shows a single component only, so the 1550 line has 
no significant effect on C I 1656 at z = 1.77637. 

• CI 1560 is covered by both the Keck and VLT spectra. It has Al III 1862 nearby, 
at z = 1.32535-1.32884, and the corresponding Al III 1854 is blended with CIV 1548 
at z = 1.78586-1.78719. The corresponding CIV 1550 shows that Al III blending of 
CI 1560 is not significant. CI* 1560 is just shortward of Fel 2484 at z = 0.74461, 
from a well-known system which shows multi-component Mgll. Fel 2523 & 2719 were 
fitted simultaneously at this redshift, showing a single velocity component which does 
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not affect the CI* lines. 

• CI 1328 and lines at shorter wavelength have coverage only from UVES at the VLT. 
The fit to C I 1328 includes a broad weak Lya at z — 2.03594, but this changes the 
effective continuum for the CI* lines mainly. It is much broader than the CI* lines 
measured. CI** 1329 is in the wing of Lya at z = 2.03696, again with small effect on 
the line parameters. 

• CI and CI* 1280 are blended with strong Lya lines at z = 1.92398 -1.92471. Only 
the ground state CI 1280 at z — 1.77637 is well constrained. 

• CI 1277 is blended with weak Lya at z = 1.91745, 1.91782 and 1.91837, affecting 
mainly the CI* and CI** lines with redshifts z > 1.7765 and rest wavelengths 
> 1277.5 A. 



A further complication is that the isotope shifts for 13 CI relative to 12 CI lines at 1656, 
1560 and 1328A are, respectively, 0.65, -3.30 and -1.65 kms -1 . These will be particularly 
important where the lines are intrinsically narrow. For this reason 13 CI was included as 
a separate spe cies, with wavelengths and oscillator strengths, for those three lines from 



Morton! (120031 ). For the other two transitions used to fit the CI, at 1280 and 1277 A, the 
spectral S/N is lower so the absence of a 13 CI component makes little difference to the final 
result. We verified this by noting that the removal of 13 CI 1328 (where the S/N is higher 
than for the 1277 or 1280 lines) from the line list had an insignificant effect on the results. 

The results of the profile fitting are shown in Fig [2J A satisfactory fit to the data 
required the presence of at least three CI velocity components, at velocities of —5.4, 11.3 
and 23.8 km s -1 relative to an arbitrarily chosen reference redshift of z = 1.77642. We 
henceforth designate these as components Nl, N2, and N3 respectively. Details of the fitted 
parameters for the C I lines, and any others required to fit the data in the regions of those 
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lines, are given in Table [U 13 CI was marginally detected in only the narrow component, 
N2. 

The CI Doppler parameter, bci = 5.08 ± 0.24 km s _1 , for component Nl corresponds, 
if the broadening is purely thermal, to a temperature of over 16 000 K, so it is likely that 
bulk motions are the major contributor. We cannot determine if it is a single component 
with some bulk motion across or within it, or if the apparent width arises because more 
than one component is present with separations less than the instrument resolution. The 
populations of the three C I levels are not consistent with a single excitation temperature - 
the best fit is T ex = 11.24 ± 0.34 K, but the observed CI** column density is then too high 
relative to the fit by A log iV = 0.65, i.e. over 4.5a. 

The unresolved CI in component N2 with bci = 0.55 ± 0.13 km s _1 has, for thermal 
line broadening, a kinetic temperature T = 220 K, with a 2a upper limit T < 480 K. This is 
much lower than the 21cm spin temperature for the whole complex, and much higher than 
the C I*/ C I excitation tem perature (which, for this component, is 6.90^gg2 K, in agreement 



with 



Songaila et al. 



( 119941 ) to within the errors) 



Component N3 is required only for the strongest two CI lines (at 1560 and 1656 A). It 
is very broad, and all we can say from its presence is that there are likely to be a number 
of weak C I components present over a velocity range ~ 50 km s _1 centered on the redshift 
given. They do not strongly affect the derived parameters for the two narrower systems. 



5. The narrow CI component 

5.1. Doppler width determination and its reliability 

The C I lines in component N2 are unresolved, and their profiles in the spectral data are 
close to the instrument profile, so we should consider further whether or not these Doppler 



13 



parameters and their error estimates are reliable. If all the lines are unsaturated then the 
Doppler parameters are not well constrained, as demonstrated by 



Narayanan et al 



(12006( 1 . 



However, it is possible to infer the Doppler parameters for the unblended C I lines to quite 
good precision provided that not all of the lines are on the linear part of the abso r ption 



Jorgenson et al. 



(120091). 



line curve o f grow th, as was demonstrated in a similar context by 
Stromgrenl ( 119481 ) describes the basic technique as applied to line doublets, and this 
generalizes to m ultiple lines, an d multiple species. T here are many other e xamples of its use 



in this way, e.g. 



Morton! ( 119751 ) 



McCandlissI ( 120031 ) and 



Cui et al 



( 120051 ). For unresolved 



lines the minimum x 2 profile fitting technique matches the equivalent widths of the lines, 
and so if there are a number of saturated as well as unsaturated lines present the Doppler 
parameter b may be determined even where there are possible blends with other components 
or ions from different redshift systems. 

To illustrate how it applies for the narrow CI component discussed here, and to 
highlight the circumstances under which we may derive reliable Doppler parameters and 
column densities from unresolved absorption lines, we have generated a new continuum 
which contains all absorption lines fitted apart from those of component N2 ground state 
CI. We then determined equivalent widths, with approximate error estimates, for the 
component N2 CI 1656, 1560, 1328, 1280 and 1277 A transitions, and show these in Fig [3] 
as error ranges on a curve of growth where the C I column density is the best fit value from 
Table HJ 



It is clear from Fig. [3] that the measured transitions with highest and the lowest 
oscillator strengths, i.e. 1656 and 1280, are vital for constraining the Doppler parameter. 
The free parameter along the x-axis is the CI column density, and had the strongest line, 
C 1 1656 not been measurable (e.g. through blending) then the remaining points adequately 
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fit the curves for any larger Doppler parameter - the CI column density could then be a 
factor of two lower, and the points on the linear part of the curve of growth, and so the 
Doppler parameter constrained only by the much larger instrument resolution value. 

If the weakest line used, C 1 1280, had not been measured, then there are much greater 
uncertainties in the column density, which could then be well over a factor of ten higher with 
the consequence that the Doppler parameter would be somewhat lower. This would reduce 
the temperature upper limit, and make any relative abundance analysis very uncertain. 

These may be obvious points to make, but when using automated profile fitting 
procedures one can forget about their inadequacies, or ignore the large error estimates. It 
is important to check that the constraints one obtains are based on the data, and not e.g. 
some approximation made in estimating the errors in the parameters. 

The VPFIT program provides error estimates from the diagonal terms in the covariance 
matrix, and these are generally reliable for systems in which the lines are resolved. For 
systems where the line widths are considerably less than that of the instrument profile the 
reliability of the estimates for the parameters and their errors has been little investigated. 
We have done so here by generating 84 simulated spectra with the same S/N as the original 
data using the line parameters given in Table [TJ For each of these Voigt profiles were 
fitted using the same fitting regions as in the original data, and the resultant parameters 
compared with those used to generate the spectra. 

From these trials the distribution of both the log column density and Doppler 
parameters alone shows that there is no significant difference between the input and 
mean output values. For the column density the mean and standard deviation are 
logiV(Cl) = 13.077 and s = 0.104, compared with an input value of logJV(Cl) = 13.064 
and VPFIT error estimate of 0.134. For the Doppler parameter input b = 0.553, with 
VPFIT error 0.125, the mean estimate from the 84 trials is b = 0.561, with a standard 
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deviation of 0.082. Therefore the difference between mean values from the simulated 
spectra differ from the input values by about the estimated error for those means (for 
the log column densities the difference is 0.013 and error estimate 0.011; for the Doppler 
parameters the corresponding quantities are 0.009 and 0.008). 

From the same trials we find that the carbon isotope ratio and its error estimate for 
this component log (iV( 12 Cl)/iV( 13 Cl)) = 1.28 ± 0.31 are reliable. A la lower limit for the 
12 C/ 13 C ratio is 5. 

Further checks were undertaken to verify that there are no additional uncertainties in 
the analysis of the data: 



The stopping criterion for the iterations involved in the fitting process was chosen 
to be when a change in y 2 < 1.0 x 10~ 5 . For any search where convergence to the 
solution is slow this could result in e.g. the Doppler parameter being biased towards 
the initial guess relative to the true value. Tests showed this concern to be unfounded 
for the stopping criterion adopted. 

Tests were run to verify that the adopted instrument resolution did not significantly 
affect the results. It is difficult to see how the resolution in each case could be worse 
than the slit-limited values which were adopted for the fits described above, but in 
conditions of good seeing the values of the FWHM appropriate for the spectra could 
be lower. We attempted to estimate the instrumental resolution by minimizing \ 2 
for profile fits to the Fell 2344, 2374, 2382, 2586 and 2600 A lines at z = 1.328, and 
found somewhat broad minima with FWHM~ 5.8 km s _1 for both the HIRES and 
UVES data. Adopting this value yields log N (CI) = 13.07 ±0.18 with b = 0.47 ±0.18 
for system N2. 



Jenkins fc Tripp! (120011 ) have estimated CI oscillator strengths from interstellar 
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medium absorption in early-type stars using HST STIS data. Jenkins (private 



commun ication) has a. 
given by 



fits with the 



so suggested that the CI 1560 oscillator strength, which is 



Morton! (2003) as fi k = 0.07 74, could be as high as f ik = 0.1316. Voigt profile 



Jenkins fc Tripp 



from those obtained using the 



20011) f - value s give results which differ by only a little 
Morton! (120031 ) values, with b = 0.53 ± 0.13 km s -1 and 



logiV(Cl) = 12.96 ±0.17, though the x 2 statistic for the fit is now somewhat too high. 



From all these trials we conclude that the parameter estimates for system N2 from 
VPFIT, for this mixture of saturated and unsaturated lines, are reliable, though the VPFIT 
error estimates are too high by a factor of ~ 1.3 - 1.5. In particular, the best estimate for 
the CI Doppler parameter from the trials is bci = 0.553 ± 0.082. This corresponds to a 
temperature of T = 220 K if the broadening is thermal, and in any case the 2a upper limit 
to the kinetic temperature for this component is 480 K. Consequently the gas in component 
N2 is very likely to be a CNM. 



5.2. Physical conditions derived from C I 

The closeness of the CI excitation temperature 6.901qq2 K to the microwave 
background temperature at the system redshift, Tcmb = 2.725(1 + z a bs) = 7.57 K, allows 
us to place an approximate upper limit to the density in the region, assuming that its 
temperature is significantly higher, since the density must be low enough that collisional 
processes are unimportant. We find a hydrogen number density nn ^ 3 cm -3 . Under 
these circumstances, if most of the elements are neutral, we can put a lower limit on 
the size of the cloud by estimating the amount of neutral hydrogen associated with the 
cloud as follows: Assuming solar relative abundances, [C/H] = —3.61, and therefore, 
logn(H) cloud = logn(CI) c i ou d + 3.61. Summing over the ground and excited states gives 
the total logn(CI) c ; OU( 2 = 13.10 cm -2 , therefore N{BI) doud = 10 16n cm -2 . The cloud must 
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therefore be greater than i = iV(HI)/n(HI) 1.7 x 10 16 cm = 0.006 pc. This lower limit 
is sufficiently small that we would, without evidence from H 2 , be concerned that the CI 
absorber only partially covers the background source. We believe this is not the case, for 
the reason given at the end of the next section (jSJ). 



6. Molecular hydrogen 

6.1. Results from a single component fit 

Further evidence for CNM gas in DLA1331+170 stems from the detection of H 2 
absorption lines at redshift z = 1.776553 in the Lyman and W erner bands with rest 



20051 ). These authors found the 



wavelengths from the Lyman limit to ~ 1120 A ( Cui et all 
rotational level populations between J = and J = 5 to be consistent with a Boltzmann 
distribution characterized by a single excitation temperature, T ex =152±10 K. This differs 
from the Galaxy ISM where low values of T ex apply for J=0 to J=2, or J = 1 to J = 3, 
while higher values of T ex apply if higher values of J are considered. The standard 
interpretation of the ISM results is that for ortho- (odd J) and para- (even J) H2 separately 
the rotational states with the lower values of J are collisionally populated which drives 
T ex toward the kinetic temperature, T, while rotational states with higher values of J are 
populated both by UV pumping and by the formation processes on grain surfaces which 
leave H 2 in the J=4 state (Spitzer & Zweibel 1974). The radiative processes are dominant 
at high J values because radiative self-shielding is less important. As a result, the constant 
value of T ex in the case of DLA1331 + 170 implies that T ex =T and that collisional processes 
dominate radiative processes for all of the rotational J levels J=0 to ,1—5. 

If collisional processes dominate for J = 4 and 5, then we can place a lower limit on the 
gas density. It must then exceed the critical density n CT it=^-j,j-2l I( 1j,j-2, where Aj^_ 2 and 
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qj,j-2 are the rates of de-excitation due to spontaneous photon decay and collisions through 
transitions between the J and J — 2 states. For the J=4— >2 transition, n crit must exceed 
~ 5 x 10 3 cm -3 when T=150 K. The gas giving rise to H 2 absorption has some fraction of 
the total iV(H I) column density associated with it, so for this region iV(H I) < 1.5xl0 21 
cm" 2 . The corresponding length-scale of a cloud with this volume density is d < 0.1 pc, 
so much less than the size of the optical continuum source of the background quasar, 
which probably exceeds 0.5 pc. Therefore, unless the gas is in a thin sheet perpendicular 
to the sightline, the smaller H2 absorbing cloud cannot cover the continuum source, which 
it must do in order to produce the saturated and damped absorption lines that arise in 
DLA1331+170. 



For these reasons 



Cui et al. 



( 120051 ) assume that radiative rather than collisional 
processes govern the population of the rotational states up to J > 4. But in that case the 
similarity between the T ex computed for states with J > 4 and J < 4 (where T ex =T) would 



Cui et al 



fl2005f ) find that at A 



be a coincidence. To obtain the required level populations, 
~ 1000 A, J v = 6.7xl0~ 23 ergs cm -2 s _1 Hz -1 sr -1 , which they claim is a reasonable value 
for FUV background radiation at z = 1.77, since if one assumes a frequency dependence of 



v 



-0.5 



J v is in accord with the value of 7.6x10 23 ergs cm 2 s 1 Hz 1 sr 1 at the Lyman 



limit inferred from the proximity effect in the Lya forest. However, the problem with this 
argument is that a power-law extrapolation across the Lyman limit frequency, while valid 
for the quasar contribution to the FUV background, ignores the dominant con t ributi on 
of Lyman Break galaxies at FUV frequencies. According to lHaardt fc Madaul (120011 ). 



J„ ~3xl0 20 ergs cm 2 s 1 Hz 1 sr 1 at A ~ 1 000 A, which is a f actor of 500 higher than 



the 



Cui et al. 



(120051 ) estimate. The effect of the lHaardt fc Madaul (120011 ) background on the 



level populations is shown in Fig. HJ which plots the quantity Nj/gj versus (Ej — E )/k, 
where Nj and gj are the column densities and degeneracies of the J th rotational state, 
and Ej is the corresponding energy eigenvalue. We show the odd and even J lower level 
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populations at the same excitation temperature since that is what the data indicates, but 
have not considered processes linking the two. 

As a result our analysis of the H2 absorption lines arising in DLA1331+170 has raised 
the following dilemma: If the J=0 — > 5 rotational states are populated by collisions, then 
the scale length of the H 2 absorbing gas would be too small. On the other hand if the 
rotational states with J > 4 are radiatively populated by plausible FUV radiation fields, 
the predicted values of Nj/gj would considerably exceed the values inferred from a single 
H 2 cloud model fit. 



6.2. Interpreting the H 2 data: Multiple component molecular hydrogen 

The presence of velocity components Nl, N2 and N3 hints at a possible solution to this 
dilemma; namely, the bulk of the H 2 gas resides in the narrow- lined component N2, with 
the broader components Nl and N3 containing a smaller amount. Consider the evidence. 
First, our best estimate for the kinetic temperature of component N2, T ~ 200 K (from 
CI if the line broadening is th ermal) , is c onsist ent with the H 2 excitation temperature, 



T ex = 152 ± 10 K obtained by 



Cuietal 



(120051 ) . Second, the redshift difference between 



component N2 and H 2 absorption, 3.0 km s _1 , is less than the difference between H 2 
and either of the other two C I components. Although the Doppler parameter for the H 2 
absorption complex, b = 13.9 ± 0.5 km s _1 is considerably greater than that of component 
N2, we suggest that while component N2 is embedded in the more turbulent medium 
that dominates the velocity structures of the H 2 absorption lines with higher J values, 
it none-the-less contains the bulk of the H 2 gas. Consequently, while component N2 is 
responsible for most of the damping-wing absorption produced for the J=0 and 1 lines, its 
low Doppler velocities imply that the absorption from states with higher J values will be 
dominated by absorption arising from the more turbulent velocity structures Nl and N3. 
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Stated differently, the values of Nj for J > 3 may have been under-estimated by the failure 
to identify the low equivalent widths generated by component N2. As a result, the Nj/gj 
ver sus (Ej — En) Ik curve m ay in fact be consistent with FUV radiation intensities due to 



the lHaardt fc Madaul (j200ll ) backgrounds and possibly low rates of local star formation. 
Finally, it is well-known that because C I and H2 are photoionized and photo-dissociated by 
photons of similar energy, they are likely to be co-spatial and hence show similar velocity 
structure. 

To investiga te this possibility further we have performed our own analysis of the 



Ouietal 



(120051 ) STIS data, assuming a component structure suggested by the CI results. 
We have taken three components for H2, two at the redshifts corresponding to the CI lines 
in components Nl and N2. Since the redshift of the third CI component, N3, is not well 
constrained the third H2 component redshift (denoted N3', z = 1.7767176 ± 0.0000067) 
was determined by the profile fit to the molecular hydrogen lines. The Doppler parameters 
were constrained so they are physically consistent with those of the corresponding CI 
i.e. bci < &h 2 < v6£>ci> with the limits corresponding to purely turbulent and purely 
thermal broadening. For the narrow component the adopted Doppler parameter was 
&h 2 = 1-0 km s _1 . The parameters for the profile fits obtained with these constraints are 
given in Table EJ which also lists the transitions used in the fits. In this case the low J 
absorption is dominated by component Nl, while the higher J is dominated more by the 
other two velocity components. This is illustrated in Fig. [5] using selected transitions for 
each J level. 

Our three-component fit to the H 2 lines is by no means unique, but it does serve to 
provide a more physically consistent picture of the complex, at least within the rather large 
errors. The single- component fit giving an apparent single excitation temperature for all 
excitation levels arises simply because different components dominate different J-values, 
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and is unlikely to reflect the true physical conditions in the individual components 
The molecular hydrogen results also suggest that here, as in the case of 



Jorgenson et al. 



( 120091 ). the narrow components inferred for CI are unlikely to be an artefact caused by a 
small covering factor for the absorbing gas, since the corresponding saturated H2 have zero 
residual intensities. 



7. Ionized heavy elements 
7.1. Singly ionized species 



While DLA1331 + 170 presents clear evidence for the presence of cold gas, it is difficult 
to determine how much of the neutral hydrogen is associated with it, and so what the 
overall cold gas fraction of the DLA may be. In an attempt to clarify the overall structure 
of the DLA we analyze the non-neutral heavy elements. 

Some of the HI will be associated with singly ionized heavy elements e.g. CII, Mgll, 
Sill, Fell etc., so that at least part of the component structure of HI will be traced by 
these ions. For the z = 1.77642 complex discussed here, the dominant components of 
all the available lines of CII and Mgll are either strongly saturated or too weak to be 
measurable, so they provide little useful information. For this complex Sill 1808, SII 1250, 
1259, Mnll 2576, Fell 1608, 2249, 2374 and Nill 1370, 1709, 1741, 1751 have the right 
combination of oscillator strength and column density to be useful, and the saturated lines 
Sill 1260, 1304, 1526 were also used to provide some additional constraints. Some other 
potentially useful lines were found to be blended with strong lines from other redshift 
systems, so were omitted from the analysis. All the elements used, with the exception 
of sulphur, have ionization potentials for I-II in the range 7.4 - 8.2 eV, and II - III from 
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15.6 - 18.2 eV, so they should arise predominantly in the same regions. For sulphur the 
corresponding potentials are 10.4 and 23.3 eV, so there could be some differences but we 
failed to find any. For AlII there is only a single transition available, at 1670A, and it was so 
saturated that the column densities of the dominant components could not be determined. 

Since we are interested in physical entities, we have fitted the multiple ions to common 
redshifts, and constrained the Doppler parameters assuming that bulk motions have a 
Gaussian distribution, and so the bulk and thermal motions for each ion add in quadrature. 



Then, for each component, b = \fb\ nTh + 2kT/m, where 6 tU rb is the turbulent (bulk) 
component, k Boltzmann's constant, m the ion mass and T the temperature. In terms of 
the VPFIT program used, this involves fitting 6 tur b and T for at least two ions of different 
mass at the same redshift simultaneously, with the constraints that both variables are 
non-negative. The error estimates from the program then apply to 6 t urb and T, not to the 
Doppler parameters for the individual ions. 

From the fitting analysis 9 separate velocity components are needed, labelled SI - S9 
in increasing redshift order. Details of these are given in Table [31 along with parameters 
for other systems with lines in the fitting regions used. All the lines fitted were available 
in the UVES spectrum, and for the HIRES data Fell 2249, Sill 1808, NiII1751, Nill 1709, 
Fell 1608 and part of Sill 1526 were included as well. The fitted line profiles and component 
structure are shown against the UVES data in Figs [6] and [71 An additional component at 
z = 1.7748722 was included, since Sill 1260 from that system is blended with SII 1259 from 
the main complex. Its structure was determined using CII 1334, Mgll 2796, 2803, Sill 1193, 
1304, 1526. Fell is not seen in this component, with an upper limit logiV(Fell) < 11.5 
(cm~ 2 , 2a). 



The componen t 



struc ture given in Table [3] is somewhat different from that given by 



Prochaska &: Wolfe! (119971 ). who illustrate a fitted profile to the Sill 1808 line. A large 
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part of the difference probably comes about because we have used the additional UVES 
spectrum and fitted many transitions where they used just one. 

For determining which ions provide the strongest constraints it can be useful to have 
Doppler parameter error estimates for individual ions rather than those for the variables 
&turb and T. Indicative values are given Table El Here the error estimates given for the 
individual Doppler parameters for each ion were computed by assuming that the best-fit 
6-value is the tied one, but computing the errors as if the tied constraint were not applied. 
In some cases the error exceeded the Doppler parameter, so there is no useful constraint. 
These are indicated by a dash in the table. 

As can be seen from Table [3j for the component structure inferred here, bulk motions 
dominate in several cases. These are S3, S4, S5, S6, S8 and S9, where the Doppler 
parameters have similar values for all ions. For SI and S7, however, thermal broadening is 
more important, though the temperature estimate for S7, at 95 000 K, is somewhat higher 
than we might expect. There are several components across the velocity range, and so, 
partly because most of the systems are now blended and partly because the mass range over 
which reliable Doppler parameters can be obtained is at most from 28 Si to 56 Fe i.e. a factor 
of less than 2, the error estimates for both b tUT \> and T are rather large. For component S7, 
the la error estimate is almost 70 000 K, so within the errors the temperature could be 
only a few xlO 4 K. 

Most temperature errors are smaller, and in two cases (S4 & S5) the 2a upper limit 
indicates temperatures of ~ 2 x 10 4 K or less. In general the temperatures are at least 
consistent with the notion that the material is in the WNM with temperatures of ~ 1 - 
3 x 10 4 K. However, the errors are large, and so CNM gas is an equally viable interpretation. 
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7.2. More highly ionized heavy elements 



We may explore the structure of more highly ionized species, though because they have 
a range of ionization potentials, we cannot necessarily treat them as coming predominantly 
from the same physical regions. The motivation for doing this is to check whether or not 
comparison of the column densities of close ions, such as Al II and Al III, provides a reliable 
indicator of the ionization state of the gas. 

The doubly ionized species detected in our observed range are Si III 1206 and Al III 1854, 
1862. Si III has only a single line which is saturated, in the Lya forest and close to the 
strong damped Lya associated with the system, so it provides little useful information. 
Fe III 1122 is potentially observable, but in a region in the Lya forest (< 3120A) where the 
S/N in the data is poor. Consequently we have used only the Al III doublet to investigate 
the component structure here. The results are given in Table HJ and the component 
structure shown against other ions in Fig. |HJ 

Similarly we can investigate the structure of triply ionized species, such as Si IV and 
CIV through their doublet absorption lines. In this case we choose Si IV only, since its 
ionization potential (45.1eV) is less than that of Hell (54.4eV), whereas that of CIV 



(64.5eV) is above the He 
in the same regions. The 



I value, and therefore the two species may not be 



p redom inantly 



Prochaska fe Wolfe! f)1999h and Wolfe fe Prochaskal (120001 ) studies 



also show that while there is a general similarity in the Si IV and CIV velocity structure, 
the correspondence is not perfect. The results for the Voigt profile fits to the Si IV 1393, 
1402 doublet are given in Table HI 
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8. Comparing the various components 



8.1. Comparing profiles 



Figure [8] shows the HI 21cm profile with representative lines from neutral, singly 
and more highly ionizied heavy elements. This is quite instructive, since there are some 
features which may be common to more than one ionization level and others which seem to 
correspond less well. 

The HI 21cm profile stands out as having little in common with any of the others. 
While we have no direct measure of what the velocity distribution of the damped Lya 
absorption is from the UVES spectrum, if the heavy element abundances are not too 
disparate between velocity components it would be surprising if it did not roughly follow 
that of the the higher column density heavy element neutral and singly ionized species. 
Therefore we expect it to follow S4 - S8 reasonably well, with some contribution from Nl 
and N2. 

The HI 21cm line is broad and, at the S/N available, smooth, and is consistent with 
the absorption arising in a single cloud at a temperature of T = 15000 K. The difference 
in shape from that expected from the heavy elements, and the redshift of the maximum 
absorption from any of the components, make it difficult to escape the conclusion that it 



and the Lya absorption have littl e to do with each ot 



ler. This reinforces the notion, raised 



in the context of other objects by 



Curran et al. 



( 120071 ) and 



Kanekar et al 



( 120091 ). that the 



21cm and Lya absorption lines do not arise in the same region(s), possibly because the 
radio source is extended and so it and the optical source are not coincident. Consequently 
the derivation of a spin temperature by using the radio 21 cm and the optical Lya lines may 
well be inappropriate. 
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The various ionization stages of the heavy element lines show considerable overlap, 
especially for the lower ionization states. However, from Voigt profile fitting each ionization 
stage independently, it is not always clear for adjacent ionization levels what velocity 
structure they have in common. This is most readily seen from Fig. |8l with details of 
redshift and Doppler parameters in Tables CD, [3] and HJ Comparative details extracted from 
these tables are given, in redshift order, in Table EJ In this table we have chosen Si II as 
the representative of the singly ionized components, since the silicon and aluminium atomic 
masses are close in value and so the Doppler parameters of any cospatial components should 
be nearly the same. 

In Table [5] possible goupings of systems are indicated as those between horizontal 
lines, and for those where the redshifts and Doppler parameters are compatible right hand 
brackets. Others within groupings may have significant components in common, especially 
where the Doppler parameters are significantly different. And, of course, the component 
structure inferred by assuming the minimum number of components consistent with the 
data may not reflect reality, and in any case does not necessarily lead to unique results 
(see Kirkman & Tytler, 1997, for a discussion applied to the Lya forest). However it is not 
immediately clear what assumptions we should make to attempt to disentangle individual 
physical components. 

One possibility is to assume the absorption for each component comes from a single 
region or cloud. If material with different ionization levels is cospatial, then the redshifts 
should be the same and the Doppler paramete rs should reflect t h e sam e bulk flows and 
temperatures. A variant of this was adopted by 



Milutinovic et al. 



(120101 ). who assumed that 



the same redshifts and Doppler parameters applied for neutral to doubly-ionized species in 
a z a bs ~ 2 sub-DLA. While this approach might be appropriate for high 6-values, we doubt 
its validity for narrow lines where thermal broadening may be important since a y^T/m 
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term which depends on the atomic mass m will make a significant contribution. Even 
allowing for this, there may well be temperature structure linked with ionization structure 
within a cloud which receives ionizing radiation from an external source. This may provide 
an explanation for the significant difference in Doppler parameters between components N2 
and S5, which share a common redshift so could well be parts of the same structure. If 
there is any AlIII associated with this component, it can not be separated from the much 
broader D4. 

Similarly, it is also tempting to suggest that Nl and S4 are parts of the same 
cloud, though the Doppler parameters are rather large and difficult to reconcile with a 
turbulent+thermal model. Also, here the component velocities relative to the reference 
redshift, Av, are significantly different. However, the velocity difference (1.1 ± 0.2 km s _1 ) 
i s small, and consistent with that which might be expected from a photdissociation region 



(jHollenbach fc Tielens 



19991 ). implying that the physical association picture may still be 
correct. So in this case a thermal and bulk flow separation within the entity could provide 
an explanation. Once again, any associated Al III is not discernible, but could be hidden in 
components D2 and D3. 

Despite the above considerations, we further explored possible constraints by forcing 
different ionization levels to have the same redshifts, but allowing the Doppler parameters 
to differ between them. In the end this is not useful because the blending of components 
leads to ill-constrained solutions, and the column densities in for the components forced into 
the blends have large errors. For example, because of the small but significant difference 
in the redshift centroids between Nl and S4, the singly ionized component S4 had to be 
divided into two, one corresponding to the CI component Nl and another, S4', still close to 
the initial S4. The column density errors in the new part of S4 corresponding to Nl were 
high (generally > 0.2 dex from VPFIT, so likely to be underestimates), and provided little 
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new information. Similarly the errors in the Al III column density in the component forced 
to the redshift of S4' were high, formally a factor of 1000, and so not useful. Some attempts 
on subsets where the Doppler parameters were also constrained yielded similar results - the 
column density errors were again too large for the values to be useful. 

While it is clear that some Al III is associated with the singly ionized components, 
some is also predominantly associated with Si IV. A similar association of some of the 



AlIII with higher ionization components was noted in GB1759 + 7539 by 



Prochaska et al 



( 120021 ). A consequence of this is that associating the total Al III column density with the 
AlII (or its proxy, Sill) will lead to an over-estimate of the amount of ionized material 
associated directly with the HI regions. Unfortunately for DLA1331 + 170 we were unable 
to associate velocity components across ionization states to obtain good relative column 
densities. Similar ambiguities arise in a detailed study of a nearby multi-phase Galactic 



interstellar cloud (INehme et al 



20081 ) for similar reasons - there are too many indistinct 



components in velocity space. 



8.2. Cold gas 



Wolfe. Gawiser & Prochaska! (120031 ) have shown that the bulk of the medium showing 
C II* absorption must be cold in general. Since we have another indicator of cold gas in C I 
which is largely distinct from that traced by the singly ionized heavy elements, we could, 
at least in principle, check this directly. In DLA1331 + 170 the CII* 1335 transition is 
blended with a Lya forest line at z = 2.0506465, so the weaker CII* 1335 lines are partly 
masked out. There are sharper CII* 1335 components within the Lya profile, and they are 
compatible with both the N1-N3 and S2-S9 structure. The C II* 1335 /Lya profile is shown 
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in iJorgenson. Wolfe fc Prochaskal (120101 ) . 



Fitting models based upon the velocity structure of the neutral and singly ionized lines, 



along with the Lya forest line, results in each c ase in summed column densit i es for ClI* 

Wolfe. Prochaska fc Gawise] ( 2003 ). They 



which are smaller than that quoted in Table 1 in 
quote logiV(Cir) = 14.05 ±0.05, while we find that, for the sum log JV(CII 
(neutral CI model) or 13.57 ± 0.55 (singly ionized model). 



13.59±0.09 



The relationship between the C I and is a 



temperature of 152 ± 10 K given by ( jCui et al- 



so no t completely clear. The H 2 excitation 



20051 ) is consistent with the kinetic 2a 



temperature limit of < 480 K we find for the C I associated with component B. However 
their estimate is based on a single-component model, while the CI which we expect should 
be co-spatial, shows at least three velocity components. We have shown in section that a 
three-component fit for the H 2 is compatible with the C I and does, within the rather large 
errors, allow the H 2 upper level populations to be consistent with the FUV background 
radiation. The temperature results are summarized in Table EJ The para- and ortho- H 2 
excitation temperatures agree at about the combined la level, though there may be a small 
radiative term enhancing the J = 3 level. 

Because the isotope shift for 13 C relative to 12 C is potentially measurable, 13 C was 
included in in the analysis, though it made little difference to the final results. Our lower 
limit of 12 C/ 13 C> 5 ( 2a) for component N2 i s considerably smaller than the lower limit 



12 C/ 13 C> 80 given by 



Levshakov et al. 



(120061 ) from their analysis of the C I 1560 and 1656 



ground state and fine structure transitions at redshift z = 1.150789 towards HE0515 — 4414. 
We suspect that the errors they give may be optimistic, since our reanalysis of the archive 
data for that object indicates that 12 C/ 13 C> 22, with the main constraint coming from 
the C I* fine structure lines. However, we are still some way from being able to use carbon 
isotope abundance limits to investigate chemical enrichment mechanisms at high redshifts. 
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Prantzos. Aubert &: Audouzd (119961 ) have shown that for timescales < 5 x 10 9 years this 
ratio is > 100. For very low metallicity stars, which may b e more appropriate for the 



DLA1331 + 170 case, that ratio can be even higher, > 1000 (ILau. Stancliffe fc Tout 



20071 ) 



Given the upper limit on the te mperature, and t h e upp er limit on the density of 



nn ~ 3 cm 3 we could consider, as in 



Jorgenson et al. 



(120091 ). whether or not the narrow 



component cloud could be gravitationally confined. However, it is not clear how much 
of the HI is associated with the cold component here, so we are not able to draw any 
firm conclusions. If we guess that, say, 10% of the HI seen is associated with the CI, 
corresponding to logA^(HI) = 20.17 (cm -2 ) then we find that the cloud is unlikely to be 
gravitationally bound. 



9. Conclusions 



We have presented here a second PL A with a hidden, c old co mponent revealed through 



narrow CI lines. The first, described by iJorgenson et al. 



(120091 ). used curve-of-growth 



techniques on an isolated CI component. Here we have shown that, even if there is a small 
amount of blending, it is possible to determine line widths, and temperature limits (in this 
case < 480 K), well below spectral resolution if the spectra contain a number of lines of the 
same ion with differing oscillator strengths, provided that they span different parts of the 
classical curve of growth. This is neither new nor surprising, but it is comforting to know 
that what are usually thought of as 'profile fitting' packages don't rely exclusively on the 
line profile information, and that the Doppler parameter error estimates are reliable. 

We note that while the spin, excitation, and kinetic temperatures within the absorbing 
system differ considerably those temperatures that should be physically consistent indeed 
are, as in the narrow component kinetic temperature and the lower J molecular hydrogen 
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excitation. Even given this agreement, the single velocity component analysis of the higher 
J H 2 excitation yielded level populations which are too low for the expected background 
radiation levels. This conflict is removed by decomposing the H 2 into three velocity 
components, corresponding to those found for CI. This example shows the potential 
importance of considering multi-component models to help interpret observed molecular 
hydrogen level populations. 

Attempts to determine temperatures for singly ionized components using transitions 
from ions with a range of atomic masses were not generally successful. For the strong 
components the transitions of the lightest ions, such as C II, are saturated, and so weaker 
lines from Si II and Fe II have to be used, and even for these blending adds significantly to 
the uncertainties. The weak limits we have obtained are consistent with temperatures of up 
to a few xl0 4 K. 

While common components across various close ionization levels are clearly evident, 
the detailed overall correspondence in velocity and Doppler parameter between the various 
ionization levels I - IV is not very good. In particular, there are differences between the 
close ions Sill and AlIII, with some of the AlIII associated instead with higher ionization 
material. Consequently, attempts to derive ionization corrections using ratios using 
these could well give misleading results unless there is a detailed study of corresponding 
velocity structures. Associating all the Al III with the lower ionization species is likely to 
overestimate the amount of ionized gas directly associated with the HI gas. While the 
results here apply to a DLA system, they suggest that the interpretation of abundances 
relative to hydrogen in sub-DLAs, where ionization corrections are more crucial, may 
need to be re-examined. This will not be very important for DLAs, where the average 
heavy element abundances in each system are reasonably well estimated by comparing 
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singly ionized heavy elements with HI, since the total column densities are usually well 
constrained (see e.g. the VPFIT version 9.5 program documentation, section 11.4) and the 
ionization corrections small. 

The profile and component disparity between the heavy elements which should trace 
the neutral hydrogen reasonably well, the neutral and singly ionized components, and the 
21cm absorption profile, reinforces the suggestion that the 21cm and Lya lines are only 
weakly related. This may be due to the fact that the background radio emission is more 
spatially extended than the optical source, and so the sightlines sampled are not the same. 
This leads to the conclusion that, if this is a typical example, the individual 21cm spin 
temperatures determined by comparing the 21cm and Lya lines are physically irrelevant. 
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T 




velocity (km/s) 

Fig. 1. — The HI damped Lya and 21cm absorption line profiles shown against unit continuum on the 
same velocity scale relative to a reference redshift z = 1.77642. The data are shown in black with Voigt 
profile fit to the damped Lya at z = 1.77674 shown in green. The fitted curve departs from the data at longer 
wavelengths in the base of the damped Lya line because of the presence of another (sub-DLA) absorption 
system at z = 1.78636. The 21cm line is also shown offset to the right and expanded by a factor of 10 in 
both x- and y- directions, with a grey tick mark indicating the reference redshift and a green one the fitted 
line centroid. 
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Fig. 2. — The CI lines on a velocity scale relative to a reference redshift z = 1.77642. The data are shown 
in black with Voigt profile fits to C I in dark blue, C I* in green, C I** as dark blue dot-dash (discernible only 
for the Nl component in the CI 1656 profile at ~ 5 km s _1 ), and unrelated lines at different redshifts as 
grey dots. The possible 13 C I component is shown in orange. The total fitted profiles are shown in red. The 
continuum is normalized to unity, and different zero offsets (of an integer x0.5) have been applied to separate 
the various transitions. The tick marks show the centroids of the fitted features for the CI components Nl, 
N2 and N3 (see text). The green tick marks indicate the positions for the C I* transitions in each wavelength 
region. 
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log(NfA) 

Fig. 3. — Curves of growth for Doppler parameters b = 0.25, 0.5, 1.0, 2.0 and 4.0 km s _1 
(solid lines) and for the component B best-fit value of b = 0.55 km s _1 (green). Estimated 
ranges for the equivalent widths w of the CI transitions at rest wavelengths A = 1656, 1560, 
1277, 1328 and 1280 give the to/A values ordered from right to left. For CI 1560 both the 
HIRES (red filled circles) and UVES (blue squares) equivalent width ranges are shown. The 
best estimate column density logiV = 13.06 was used to set the x-position, and the error 
range correponding to the uncertainty in logiV is shown against the 1656A line. 
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Fig. 4. — The distribution of H2 over J-levels 0-5. The abscissa is the energy of the level 
relative to J = 0, and the ordinate is t he column densit y divided by the level degeneracy. The 



black points with error bars show the 



Cui et al. 



(120051 ) values with error bars, and the green 



points are the expected values for an excitation temperature of T ex = 150 K. The red points 
show the expected distribution for a mixture of collisional excitation at this temperature 
and the Haardt-Madau background radiation at redshift z = 1.7765 for a range of densities 
as indicated, with the highest points corresponding to the lowest densities. 
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Fig. 5. — Representative H2 absorption lines from the J = to 5 levels for the three-component model 
described in the text. The data are shown in black, the overall fit in green, and the individual contributions 
from the z = 1.7763702, 1.7765246 and 1.7767176 components in red, blue and turquoise respectively. The 
zero- velocity reference is z — 1.776420, the continuum is set to 0.8 everywhere, and successive lines are offset 
in y by 1.0. The narrow component (blue) dominates the J = and 1 lines, but is only a minor contributor 
for J = 4 and 5. 
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Fig. 6. — The singly ionized complex centered on redshift z = 1.776420 showing the total fits in red. These 
include unrelated blends and the component structure for the lines in the complex. Component velocities 
and labels (see Table are marked in the corresponding colours. The continuum level is set at 0.9 and each 
line biassed up by an integer amount to separate them. Fe II 2586 was not included in the fitted regions 
because of the possibility that the weak structure seen in the blue wing may extend into the line itself, but is 
shown here as an example of the consistency of the fit obtained using other lines. The dark grey component 
shown in the S II 1259 profile is Sill 1260 at z = 1.7748722. 
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Fig. 7. — The singly ionized complex centered on redshift z = 1.776420 showing the total fits in red, and 
the component structure for the lines in the complex for the weaker lines. Component velocities and labels 
(see Table [3]) are marked in the corresponding colours. The continuum and zero level for the Fell 1608 are 
unity and zero, as shown, but for the other lines the vertical scale has been stretched by a factor of four, so 
the zero levels are four units below the continuum which has been biassed upward in each case by an integer 
amount to separate the lines. 
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Fig. 8. — The velcity structure of the HI 21cm, CI 1656, Si II 1808 Al III 1854 and Si IV 1402 absorption 
lines on a common velocity scale. The total Voigt profile fits are shown in red, and individual components 
for the various ions shown in colours corresponding to the labels given in Tables [I] [3] and 21 
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Table 1: CI and blends 



Ion 


z 


± 


b 


± 


log AT 


± 


Notes 


Fe I 


0.7446128 


0.0000006 


1.20 


0.18 


12.04 


0.02 


Fe I 2484 near C I* 1560 


Al III 


1.3253472 


0.0000014 


6.26 


0.32 


12.30 


0.02 


Al III 1862 near C I 1560 


Al III 


1.3281613 


0.0000072 


7.17 


0.86 


12.12 


0.06 




Al III 


1.3282625 


0.0000048 


6.22 


0.78 


12.21 


0.05 




Al III 


1.3285315 


0.0000014 


11.87 


0.27 


12.55 


0.01 




Al III 


1.3288447 


0.0000057 


9.65 


1.15 


11.76 


0.04 




C I 


1.7763702 


0.0000009 


5.08 


0.24 


13.08 


0.02 


Nl (-5.4 ± 0.1 km s _1 ) 


C I* 


1.7763702 




5.08 




12.59 


0.02 




C I** 


1.7763702 




5.08 




11.90 


0.14 




C I 


1.7765246 


0.0000014 


0.55 


0.13 


13.06 


0.13 


N2 (11.3 ± 0.2 km s _1 ) 


C I* 


1.7765246 




0.55 




12.05 


0.07 




C I** 


1.7765246 




0.55 




< 


11.70 




13 C I 


1.7765246 




0.53 




11.78 


0.28 




C I 


1.7766400 


0.0000467 


24.45 


6.18 


12.51 


0.11 


N3 (23.8 ± 5.1 km s" 1 ) 


C I* 


1.7766400 




24.45 




< 


12.20 




C I** 


1.7766400 




24.45 




< 


12.25 




C IV 


1.7858576 


0.0000180 


28.88 


2.38 


12.93 


0.03 


C IV 1550 near Al III 1854 


C IV 


1.7861859 


0.0000117 


10.49 


2.06 


12.43 


0.13 




C IV 


1.7864237 


0.0000116 


16.49 


1.75 


12.80 


0.04 




C IV 


1.7869053 


0.0000121 


14.83 


1.75 


12.62 


0.05 




C IV 


1.7871942 


0.0000187 


15.63 


2.70 


12.47 


0.07 




H I 


1.9174518 


0.0000518 


16.90 


9.29 


11.93 


0.23 


Lya in C I 1277 region 


H I 


1.9178188 


0.0000201 


7.20 


3.74 


11.86 


0.16 




H I 


1.9183670 


0.0000841 


10.65 


11.33 


11.86 


0.43 




?? 


1.9238783 


0.0000317 


7.78 


2.01 


12.89 


0.25 


Line in C I 1280 region 


?? 


1.9239821 


0.0000106 


6.60 


0.93 


13.23 


0.11 




H I 


1.9242860 


0.0000042 


11.48 


0.71 


13.32 


0.03 




?? 


1.9243053 


0.0000064 


1.21 


3.54 


12.86 


2.02 




?? 


1.9244645 


0.0000067 


2.21 


2.12 


12.25 


0.09 




H I 


1.9247054 


0.0000138 


9.95 


2.44 


12.18 


0.08 




C IV 


1.9661609 


0.0000093 


16.53 


1.57 


12.47 


0.04 


C IV 1550 near C I 1656 


H I 


2.0359392 


0.0000731 


37.31 


12.66 


12.14 


0.17 


Lya in C I 1328 region 


H I 


2.0369588 


0.0000458 


19.80 


7.66 


12.03 


0.24 




H I 


2.0376443 


0.0000388 


55.83 


9.36 


12.91 


0.09 





Note. — Doppler parameters b are in km s _1 , and column densities are log cm~ 2 . Error estimates are 
la, and upper limits are 2a. The C I wavelengths used were those for 12 C except where the isotope is given 
explicitly. Lya lines with b < 10 km s _1 are probably unidentified heavy elements, and are marked '??'. 
Their parameters were determined assuming that the rest wavelength and oscillator strength are the same 
as for Lya. For the C IV and Al III doublets both lines were included in the fits, and for Fe I at redshift 
z = 0.774613 the lines at 2167, 2484 and 2719 were used. Component labels are shown in bold, and the 
velocities given for the three identified components are relative to a reference redshift z = 1.77642. 
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Table 2: Multicomponent molecular hydrogen column densities 



Component Nl N2 N3' 

z 1.7763702 1.7765246 1.7767176 Lines used in fit 

6(H 2 ) (kms" 1 ) 8.7 ±0.5 1.0 - 10.0 ±0.7 



logiV(J=0) 


18.56 


±0.28 


18.99 


±0.16 


15.69 


±0.63 


1049 1077 1092 


log 7V(J=1J 


18.90 


±0.22 


19.46 


±0.08 


17.39 


±0.40 


932 1049 1051 1064 1077 
lOta 1092 


logiV(J=2) 


16.68 


±0.24 


18.44 


±0.11 


16.10 


±0.13 


927 987 1005 1016 1051 
1053 1064 1066 1079 1081 
1096 


logiV(J=3) 


15.94 


±0.14 


18.20 


±0.14 


16.07 


±0.14 


934 935 942 960 987 995 
997 1006 1017 1041 1053 
1056 1081 1096 1099 


logiV(J=4) 


14.95 


±0.06 


14.42 


±0.35 


14.81 


±0.06 


1017 1035 1047 1060 1074 
1085 1099 


logiV(J=5) 


14.46 


±0.07 


13.61 


±0.33 


14.42 


±0.07 


996 997 1006 1017 1048 
1052 1065 1079 1089 


T ex (J=0&2) (K) 


86 


+14 
-10 


177 


+30 
-22 


> 200 







Note. — Column densities are log cm~ 2 . The error estimates are lcr, and for the (log) column densities are 
likely to be underestimates if they exceed <~ 0.3. The excitation temperature for component N3' is effectively 
unconstrained - a la lower limit is given 
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Table 3: Singly ionized species and blends 



Ion 




z 


± 


b 


± 


log TV 


± 


Sys 

&turb 


Av (km s ) 
± T (10 4 K) 


± 


C II 


1 


.7748722 


0.0000014 


5.36 


(0.42) 


13.20 


0.04 


SI 


-167.1 




Mgll 


1 


.7748722 




3.76 


(0.28) 


12.11 


0.01 


0.0 


0.9 2.1 


1.1 


Sill 


1 


.7748722 




3.50 


(0.64) 


12.27 


0.03 








Sill 


1 


.7760202 


0.0000039 


6.46 


0.52 


12.70 


0.03 


S2 


-43.2 ±0.4 




Sill 


1 


.7761772 


0.0000078 


5.33 


(0.89) 


13.21 


0.10 


S3 


-26.2 ±0.9 




S II 


1 


.7761772 




5.20 


(3.50) 


13.44 


0.13 


4.3 


4.2 1.7 


5.8 


Fell 


1 


.7761772 




4.82 


(2.23) 


12.46 


0.14 








Sill 


1 


.7763599 


0.0000013 


9.53 


(0.24) 


15.07 


0.01 


S4 


-6.5 ±0.2 




S II 


1 


.7763599 




9.53 


(0.48) 


14.86 


0.01 


9.5 


0.4 0.0 


1.3 


Mnll 


1 


.7763599 




9.53 


(1.08) 


12.16 


0.03 








Fell 


1 


.7763599 




9.53 


(0.26) 


14.30 


0.02 








Nill 


1 


.7763599 




9.53 


(0.73) 


13.10 


0.03 








Sill 


1 


.7765303 


0.0000028 


5.67 


(0.99) 


14.35 


0.09 


S5 


11.9 ±0.3 




S II 


1 


.7765303 




5.67 


(1.02) 


14.32 


0.10 


5.7 


4.1 0.0 


0.6 


Fell 


1 


.7765303 




5.67 


(3.80) 


13.07 


0.27 








Nill 


1 


.7765303 




5.67 


- 


11.43 


1.01 








Sill 


1 


.7766255 


0.0000297 


22.78 


(2.34) 


14.38 


0.18 


S6 


22.2 ±3.2 




S II 


1 


.7766255 




22.77 


- 


13.73 


0.79 


22.6 


4.7 1.1 


39 


Mnll 


1 


.7766255 




22.72 


(15.29) 


12.00 


0.06 








Fell 


1 


.7766255 




22.71 


(6.51) 


13.99 


0.08 








Nill 


1 


.7766255 




22.71 


(9.38) 


12.98 


0.08 








Sill 


1. 


.7767156 


0.0000043 


8.11 


(1.53) 


14.23 


0.12 


S7 


31.9 ±0.5 




S II 


1 


.7767156 




7.66 


(2.04) 


14.03 


0.16 


3.1 


4.2 9.5 


6.6 


Fell 


1 


.7767156 




6.14 


(1.46) 


13.52 


0.11 








Nill 


1 


.7767156 




6.03 


(3.70) 


12.33 


0.15 








Sill 


1 


.7768523 


0.0000029 


8.00 


(0.55) 


14.27 


0.03 


S8 


46.7 ±0.3 




S II 


1 


.7768523 




7.97 


(2.53) 


13.99 


0.07 


7.7 


0.8 0.8 


2.0 


Mnll 


1 


.7768523 




7.86 


(6.00) 


11.51 


0.13 








Fell 


1 


.7768523 




7.86 


(0.67) 


13.87 


0.04 








Nill 


1 


.7768523 




7.85 


(1.82) 


12.72 


0.05 








Sill 


1 


.7770083 


0.0000030 


4.01 


(0.42) 


13.32 


0.04 


S9 


63.5 ±0.3 




Fell 


1 


.7770083 




3.85 


(1.48) 


12.76 


0.09 


3.7 


2.0 0.4 


2.4 


H I 


1.8757384 


0.0000243 


23.39 


0.94 


14.33 


0.09 








H I 


1.8775792 


0.0000028 


35.60 


0.48 


14.03 


0.01 








H I 


1.8760396 


0.0001075 


17.46 


6.86 


13.25 


0.53 








H I 


1.9779551 


0.0000061 


21.84 


0.88 


12.97 


0.01 








H I 


1.9799652 


0.0000161 


26.69 


2.63 


12.70 


0.03 








?? 


2.6737241 


0.0000526 


49.43 


5.17 


12.73 


0.06 









Note. - Component labels are given in bold. Error estimates are la, and those in parentheses are 
indicative Dopplcr parameter errors if the ion were treated without reference to others in the system (see 
text for details). Other details as for Table [U 
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Table 4: A1III & SilV absorption 



Ion 




z 




± 







AT 

lug 1 V 


_|_ 


bys 


A I? (km s ) 


Aim 


1.7760362 


0.0000047 


7 19 


u.yj 


1 1 QQ 


U.U4: 


Dl 


/M K _|_ n c. 

— ii .o m u.o 


Aim 


1.7763051 


0.0000172 


8.93 


1.56 


12.40 


0.11 


D2 


-12.4 ± 1.9 


Aim 


1.7764071 


0.0000098 


4.74 


1.69 


12.08 


0.21 


D3 


-1.4 ±1.1 


Aim 


1.7767114 


0.0000107 


18.85 


1.41 


12.57 


0.04 


D4 


31.5 ± 1.2 


Aim 


1.7768598 


0.0000052 


4.56 


1.67 


11.83 


0.13 


D5 


47.5 ± 0.6 


SilV 


1 


.7748504 


0. 


.0000012 


6.62 


0.16 


13.36 


0.01 


Tl 


-169.53 


SilV 


1. 


.7750657 


0. 


.0000223 


10.35 


3.50 


12.81 


0.36 


T2 


-146.27 


SilV 


1 


.7750956 


0. 


.0000062 


1.30 


1.03 


12.39 


0.14 


T3 


-143.04 


SilV 


1. 


.7752779 


0. 


.0001684 


23.70 


10.30 


12.77 


0.48 


T4 


-123.34 


SilV 


1. 


.7752931 





.0000043 


5.12 


1.13 


12.71 


0.12 


T5 


-121.71 


SilV 


1 


.7758053 


0. 


.0000093 


6.25 


1.92 


11.89 


0.08 


T6 


-66.4 ± 1.0 


SilV 


1. 


.7760357 


0. 


.0000012 


6.82 


0.22 


13.90 


0.03 


T7 


-41.5 ± 0.1 


SilV 


1 


.7762370 





.0000014 


6.52 


0.30 


13.63 


0.02 


T8 


-19.8 ± 0.2 


SilV 


1. 


.7764205 


0. 


.0000098 


8.94 


3.32 


12.47 


0.14 


T9 


0.1 ± 1.1 


SilV 


1. 


.7766329 





.0000282 


9.42 


2.71 


12.97 


0.25 


T10 


23.0 ±3.0 


SilV 


1. 


.7767771 





.0000166 


10.16 


1.95 


13.29 


0.12 


Til 


38.6 ± 1.8 


SilV 


1 


.7769879 


0. 


.0000085 


8.97 


2.69 


12.66 


0.31 


T12 


61.3 ± 0.9 


SilV 


1. 


.7773069 





.0000041 


10.41 


1.65 


13.29 


0.21 


T13 


95.8 


SilV 


1. 


.7773179 


0. 


.0000567 


29.11 


50.81 


12.96 


0.18 


T14 


96.9 


SilV 


1. 


.7774448 





.0000038 


2.17 


1.05 


12.59 


0.08 


T15 


110.6 


SilV 


1. 


.7776712 





.0000063 


5.84 


1.77 


12.52 


0.20 


T16 


135.1 


SilV 


1. 


.7778277 


0. 


.0000541 


15.81 


5.97 


12.40 


0.33 


T17 


152.0 



Note. — Error estimates are la. The velocities Av are relative to a reference redshift z = 1.77642. See 
Table Q] for further details. 
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Table 5: Velocity components compared 



Ion 




Z 




-1- 

It 







± 


log 


AT 


± 


ays 


L\U I Kill D J 




QiTV 
oil V 


i 

i 


( I OoUOo 


n 
u 


uuuuuyo 





OX. 


1 


92 


1 1 
1 1 


QO 

oy 





08 


J. D 


— 00.4 zt l.U 




O *TT 

bill 


1 


7760202 





HAAAAOA 

0000039 


6 


46 





52 


12 


70 





03 


S2 


—43.2 ± 0.4 




A 1TTT 
Allll 


i 

i 


i (OUoOZ 


n 
u 


UUUUU4 / 


7 


i f > 





95 


1 1 
1 1 


(Hi 





04 


i ) i 


— 41.0 zt U.O 






SilV 


1 


7760357 





AAAAA1 O 

0000012 


6 


o o 

82 





22 


1 Q 

13 


90 





03 


17 


A 1 C | A 1 

41.5 ± 0.1 






Sill 


1 


7761772 





0000078 


5 


33 





89 


13 


21 





10 


S3 


-26.2 ± 0.9 




SilV 


1 


7762370 





0000014 


6 


52 





30 


13 


63 





02 


T8 


-19.8 ±0.2 




Aim 


1 


7763051 





0000172 


8 


93 


1 


56 


12 


40 





11 


D2 


-12.4± 1.9 




Sill 


1 


7763599 





0000013 


9 


53 





24 


15 


07 





01 


S4 


-6.5 ±0.2 




C I 


1 


7763702 





0000009 


5 


08 





24 


13 


08 





02 


Nl 


-5.4 ±0.1 




Aim 


1 


7764071 





0000098 


4 


74 


1 


69 


12 


08 





21 


D3 


-1.4 ± 1.1 






SilV 


1 


7764205 





0000098 


8 


94 


3 


32 


12 


47 





14 


T9 


0.1 ± 1.1 






C I 


1 


7765246 





0000014 





55 





13 


13 


06 





13 


N2 


11.3 ±0.2 




Sill 


1 


7765303 





0000028 


5 


67 





99 


14 


35 





09 


S5 


11.9 ±0.3 




Sill 


1 


7766255 





0000297 


22 


78 


2 


34 


14 


38 





18 


S6 


22.2 ±3.2 


1 


SilV 


1 


7766329 





0000282 


9 


42 


2 


71 


12 


97 





25 


T10 


23.0 ±3.0 




C I 


1 


7766400 





0000467 


24 


45 


6 


18 


12 


51 





11 


N3 


23.8 ±5.1 


J 


A1III 


1 


7767114 





0000107 


18 


85 


1 


41 


12 


57 





04 


D4 


31.5 ± 1.2 




Sill 


1 


7767156 





0000043 


8 


11 


1 


53 


14 


23 





12 


S7 


31.9 ±0.5 




SilV 


1 


7767771 





0000166 


10 


16 


1 


95 


13 


29 





12 


Til 


38.6 ± 1.8 




Sill 


1 


7768523 





0000029 


8 


00 





55 


14 


27 





03 


S8 


46.7 ±0.3 






A1III 


1 


7768598 





0000052 


4 


56 


1 


67 


11 


83 





13 


D5 


47.5 ±0.6 






SilV 


1 


7769879 





0000085 


8 


97 


2 


69 


12 


66 





31 


T12 


61.3 ±0.9 




Sill 


1 


7770083 





0000030 


4 


01 





42 


13 


32 





04 


S9 


63.5 ±0.3 





Note. — Component labels are given in bold. Systems which may be associated are grouped between 
horizontal lines, and those for which the redshifts and Doppler parameters are compatible are shown linked 
with a right bracket. See Tabled] for other details. 
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Table 6: Temperature estimates 



Ion & method 


z 


± 


b 


± 


T 


± 


Notes 


H2 T ex J=o&2 


1.7765246 


0.0000014 


1.00 




177 


30 


Component N2 


H2 T ex J=i&3 


1.7765246 


0.0000014 


1.00 




225 


28 


Component N2 


CI T cx 


1.7763702 


0.0000009 


5.08 


0.24 


10.4 


0.5 


Component Nl 


CI T cx 


1.7765246 


0.0000014 


0.55 


0.13 


7.2 


0.8 


Component N2 


CI Doppler 


1.7765246 


0.0000014 


0.55 


0.08 


220 


< 480 


Component N2 



Note. - 



Error estimates are la, and upper limits are 2a. 



